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Generalized parton distributions provide information on the longitudinal and transverse distribution of partons 
in the fast moving nucleon. Furthermore, they contain information on the spin structure of the nucleon. First 
results of a lattice study of generalized parton distributions are presented. 



1. INTRODUCTION 

Understanding the internal structure of hadrons 
in terms of quarks and gluons (partons), and in 
particular how quarks and gluons provide the 
binding and spin of the nucleon, is one of the 
outstanding problems in particle physics. 

Electromagnetic form factors provide informa- 
tion on the spatial distribution of charge and 
magnetization in the nucleon, irrespective of the 
partons' momenta and independent of the reso- 
lution scale. Parton distributions measure the 
probability ^(x)! 2 of finding a parton with frac- 
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tional longitudinal momentum x in the fast mov- 
ing nucleon at a given transverse resolution 1/Q, 
while no information on the transverse distribu- 
tion of partons is provided. Generalized par- 
ton distributions (GPDs) [T] describe the co- 
herence of two different hadron wave functions 
tp'(x + £/2) ip(x — £/2), one where the parton car- 
ries fractional momentum x + £/2 and one where 
this fraction is x — £/2, from which further in- 
formation on the transverse distribution of par- 
tons can be drawn In the limit where the 
momentum transfer A on the nucleon is purely 
transverse, i.e. A = (0, Ax) and £ = r\ = 0, 
GPDs regain a probabilistic interpretation [3]. 
When Fourier transformed to impact parameter 
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• Partem density 




• Generalized parton 
distribution at r]=0 




Figure 1. Probabilistic interpretation of form fac- 
tors, parton distributions and generalized parton 
distributions at r\ = in the infinite momen- 
tum frame (resolution: 5z± = 1/Q). 



space, they measure, for example, the probability 
of finding a parton of momentum fraction x at 
the impact parameter b^. In Fig. 1 we compare 
the information provided by form factors, parton 
distributions and GPDs. In other words, GPDs 
extend the well-known Feynman parton distribu- 
tion and electromagnetic form factors of the nu- 
cleon to new kinematic dimensions. In the for- 
ward limit (A = 0), these distributions reduce to 
the Feynman parton distributions. 

Moments of GPDs are amenable to lattice cal- 
culations |5I6| . Thus, they offer promising ways 
to link phenomenological observations to first 
principle theoretical considerations. In this talk 
we will report on recent quenched results obtained 
by the QCDSF collaboration |S]. The calcula- 
tion employs non-perturbativcly improved Wil- 
son fermions and non-perturbative renormaliza- 
tion factors [7] of the operators whose matrix el- 
ements are involved. Generally, improved Wil- 
son fermions show only a very mild cut-off depen- 
dence |S]. We may therefore restrict ourselves to 
calculations done on the 16 3 32 lattice at f3 = 6.0. 

2. (GENERALIZED) FORM FACTORS 

We shall consider the GPDs H q and E„ in the 
nucleon, where q = u, d, ■ ■ ■ denotes the flavor of 
the struck quark. We will not consider the gluon 
sector here. 

The zeroth moments of H q and E q give the elec- 
tromagnetic Dirac and Pauli form factors: 

J dxE q (x,Z,A 2 ) = F2(A 2 ). 

Both form factors have been computed on the lat- 
tice |8I9| . In Fig. 2 we show a typical result for the 
Dirac form factor of the proton, F\ = + Ff. 
The solid curve is a dipole fit of the form 

Fi(A 2 ) = Fi(0)/(1 - A 2 /Ml) 2 . (2) 

In Fig. 3 the dipole mass Mi is extrapolated to 
the physical pion mass. We find good agreement 
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Figure 2. The Dirac form factor Fi(A 2 ) of the 
proton at m-n = 450 MeV, together with a dipole 
fit. 
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Figure 3. The dipole mass M\ as a function of 
together with a linear extrapolation to the 
physical pion mass (□). 



with the phenomenological value, which is ap- 
proximately equal to the physical p and u> mass. 

The first moment gives us generalized form fac- 
tors of a class of tensor operators: 

J dx x H q (x, £, A 2 ) = A\{A 2 ) + £ 2 Cf (A 2 ) , 

J' dxxE q (x,t,A 2 ) = Bl(A 2 )-ec!(A 2 ), 

(3) 

where A\(A 2 ), B q 2 (A 2 ) and C|(A 2 ) are given 
by the nucleon matrix elements of the energy- 
momentum tensor (EMT): 

(P'PI U} \P) = ^lp'\qi {l ,D v} q\p) 

= A q 2 {A 2 )u{p') 1{tl p v} u{p) 

- B q 2 (A 2 ) ^—u(pi)A a <j a ^p u} u(p) 
+ CUA 2 )—u(p')u(p)A {lM A v} . 

(4) 

Here mpj denotes the nucleon mass, p = 
| {p' + p) , A = p' — p 7 and curly brackets refer 
to symmetrization of indices and subtraction of 
traces. The EMT has twist two and spin two. 
It is assumed to be renormalized at the scale Q, 
which makes ^1(A 2 ), S|(A 2 ) and C|(A 2 ) scale 
and scheme dependent. For the classification of 
states of definite J PC contributing to (0} in the 
i-channel see |TJ]]. The so-called skewedness pa- 
rameter £ is defined by £ = — n ■ A, where n is a 
light-like vector with n ■ p — 1, and bounded by 
\g\ < 2^/A 2 /(A 2 - 4m 2 N ). In the forward limit, 
A 2 — £ = 0, the first moment of H q reduces to 
the first moment of the unpolarized parton dis- 
tribution: 

A*(0) = (x q )= f dxx( qi (x)+ qi (x)), (5) 
Jo 

where <7j( n (%) is the distribution with quark spin 
parallel (antiparallel) to the spin of the nucleon. 
Furthermore, we have ^I] 

i(Al(0) + B!(0))=J g! (6) 
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Figure 4. The generalized form factors B% 
and C% at = 450 MeV, together with the 
dipole fit and the extrapolated values at A 2 = 

(□)■ 



Figure 5. The generalized form factors A%, B% 
and C2 at = 450 MeV, together with the 
dipole fit and the extrapolated values at A 2 = 
(□)• 
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where J q is the angular momentum of the q quark, 
and 

J = J2 J i ( 7 ) 

g 

is the total angular momentum of the nucleon car- 
ried by the quarks. 

In Figs. 4 and 5 we show the generalized form 
factors Af(A 2 ), B%(A 2 ), C£(A 2 ) and Af(A 2 ), 
B%{A 2 ), C|(A 2 ) of the proton, again for = 
450 MeV. As for the nucleon form factors, the 
generalized form factors can be well described by 
the dipole ansatz 

^(A 2 )=A!(0)/(1-A 2 /Af 2 ) 2 , (8) 

and similarly for B\ and C|. The numbers refer 
to the MS scheme at the renormalization scale 
Q = 2GeV. Fits of A%(A 2 ) and ^(A 2 ) give the 
same dipole mass M within errors. We therefore 
fit our data by a common dipole mass M, only 
depending on the quark mass. For a reliable ex- 
trapolation to A 2 = it is important to cover a 
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Figure 6. The dipole mass M as a function of 
mjr, together with a linear extrapolation to the 
physical pion mass (□). 



wide enough range of A 2 . Our data do not favor 
a monopole behavior. 

In Fig. 6 we show the dipole mass M as a func- 
tion of the pion mass. Again, the numbers ap- 
pear to lie on a straight line. A linear extrap- 
olation in m,r to the physical pion mass gives 
M = 1.1(2) GeV. This value is close to the phys- 
ical masses of the ji and a 2 mesons, which sup- 
ports the hypothesis of tensor meson dominance. 
The numbers A*(0), Bf (0) and Cf (0) show little 
variation with the quark mass and are extrap- 
olated quadratically in m n to the physical pion 
mass. We obtain 



q 


Am 


Bf(0) 


cm 


u 


0.400(22) 


0.334(113) 


-0.134(81) 


d 


0.147(11) 


-0.232(77) 


-0.071(42) 



Higher moments of H q and E q may be ob- 
tained in a similar way. Knowing sufficiently 
many moments, we can reconstruct H q {x,^, A 2 ) 
and E q (x,t;, A 2 ) by inverse Mellin transform |12| . 

3. NUCLEON ANGULAR MOMENTUM 

The angular momentum decomposes, in a gauge 
invariant way, into two contributions: 

Jq = L g + S q , (9) 

where L q is the orbital angular momentum and 

Sq = l Aq ^\J dx ^ {X) " Ql {X) ^ (10) 

is the spin of the quark. We know Aq from sepa- 
rate calculations |13I14 |. so that L q can be com- 
puted from 

In Fig. 7 we show the total angular momentum 
J = J u + Jd- The dependence on the pion mass is 
rather flat, as expected JS]- The errors are due 
to the relatively large statistical errors of BV; and 
£>2 and the fact that and B\ cancel to a large 
extent. A quadratic extrapolation in to the 
physical pion mass gives 



J 


Ju 


Jd 


Su 


Sd 


0.33(7) 


0.37(6) 


-0.04(4) 


0.42(1) 


-0.12(1) 
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Figure 7. The total angular momentum J, to- 
gether with a quadratic extrapolation to the phys- 
ical pion mass (□). 



4. GENERALIZED PARTON DISTRI- 
BUTION 

If the dipole behavior (J2J , I© continues to hold 
for the higher moments as well, and if we as- 
sume that the dipole masses continue to grow in 
a Regge-like fashion, we may write 



dxx n H q (x,0,A 2 ) = «)/(l-A 2 /M 2 +1 ) 2 : 



with 



Mf 



const. + l/a' , 



(12) 



(13) 



where const, « -0.5 GeV 2 and l/a' w 1.1 GeV 2 . 
We know (a;™) from previous lattice calculations 
for < n < 3 [T7j. That is sufficient to com- 
pute H q (x, 0, A 2 ) by moans of an inverse Mcllin 
transform |12j . 

Having done so, the desired probability distri- 
bution of finding a parton of momentum fraction 
x at the impact parameter bj_ (see Fig. 1) is then 
obtained by the Fourier transform of H q (x, 0, A 2 ): 



The numbers for S q refer to our latest results |16j , 
computed from the non-perturbatively improved 
axial vector current with non-perturbative renor- 
malization factors. It turns out that the to- 
tal angular momentum J carried by the quarks 
amounts to as 70% of the spin of the (quenched) 
proton, leaving a contribution of as 30% for the 
gluons. The major contribution is given by the 
u quark, while the contribution of the d quark 
is found to be negligible, which hints at strong 
pairing effects. 

From J and S we can compute L q now. The 
total orbital angular momentum of the (valence) 
quarks turns out to be consistent with zero: 

L = L u + La = 0.03(7). (11) 

This indicates that (at virtuality Q = 2 GeV) the 
parton's transverse momentum in the (quenched) 
proton is small. 



H q (x, 0, bi) J* A ^ e ib - A -^(x, 0, A 2 ) . 

(14) 

In Fig. 8 we show a first attempt at a lattice 
calculation of H q (x, 0, bj_) for the u quark dis- 
tribution in the proton. The volume under the 
two-dimensional surface (taken over the full kinc- 
matical range < x < 1 and — oo < b < oo) gives 
us back (x u ). We see that at smaller values of x 
the distribution in b_L is rather broad (and dif- 
ferent from the distribution sketched in Fig. 1), 
while at larger values of x it becomes somewhat 
more narrow. We also notice that the peak of 
the distribution is shifted slightly towards larger 
values of x as |bj_| is increased. 

5. CONCLUSION 

We have presented a first quenched lattice 
estimate of the generalized parton distribution 
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Figure 8. The generalized parton distribution xH u (x, 0, bj_) as a function of x and the impact parameter 
b = |bj_|. The plot extends over the range 0.1 < x < 0.9 and 0.1 fm < b < 0.6 fm. 



H q (x, 0, b^_), assuming that the generalized form 
factors of the higher operators (i.e. higher than 
the tensor operator) can be described by a dipole 
ansatz as well, with dipole masses lying on a 
Regge trajectory. As a by-product we obtained 
the total and orbital angular momentum of the 
nucleon carried by quarks. A calculation of the 
missing generalized form factors is in progress. 
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